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1. Introduction 

Several recent studies have shown the advantages of using composite nucleons in finite 
density calculations PP, [2], [S!- For example, in some models the saturation mechanism 
of the system can be related to how the structure of the nucleon changes with density 0]. 
While saturating nuclear systems can also be readily obtained using point-like nucleons 
[S], if one wishes to explore the transition to quark matter it is desirable to start with a 
description of the nucleon on the quark level. Although the nucleon is a complex object, 
some simplification arises by treating it as an interacting quark-diquark state in the 
Faddeev approach. It turns out that this picture of the nucleon has considerable appeal, 
both theoretical and observational|. For example, it nicely reproduces the light baryon 
spectrum [7] , while a calculation without diquark correlations predicts an abundance of 
missing resonances |H1 . Indeed, the widespread consensus that quark pairing is favorable 
in high density quark matter [0], where perturbative QCD has some predictive power, 
lends support to the idea that diquark correlations within the nucleon are important 
too [6^. 

In this work we consider the possibility of a connection between the diquark 

correlations within the nucleon and the pairing of quarks anticipated in color 

superconducting quark matter. We work with the assumption of two fiavors. At 

zero temperature and normal nuclear matter density strange quarks certainly do not 

feature in the dynamics of the system, so this is a reasonable starting point. At 

higher densities strangeness could play a role through the formation of hyperons, kaon 

condensates and/or strange quark matter, depending on which of these phase transitions 

is favored and in what sequence they appear with respect to density (TUl HH [121 IE] §■ 

Actually, there are numerous possibilities for quark pairing and condensation in high 

density matter, especially if strangeness is introduced. Although these ideas can all be 

investigated within the framework we are using, we wish to concentrate here on the 

question of whether a unified description of diquark interactions can be achieved for 

both the nuclear matter (NM) and quark matter (QM) phases. For this purpose, we 

consider the simplest possibility that two fiavor hadronic matter goes directly to two 

fiavor QM, using the Nambu-Jona-Lasinio (NJL) model to describe both phases. In NJL 

X For a discussion on the evidence in favor of diquarks see Ref. jBj. 

§ The possibility of 3-flavor quark stars was first investigated in 1970 by Itoh |14|. 
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type models, it has been shown that the mass of the strange quark tends to inhibit the 
formation of strange QM in the density region relevant to compact stars P ^ ITH t lT7 | ll8j. 
The possibility of pseudoscalar condensation is considered in Refs.[T9 | l2( H l2 H l22]. and the 
competition between the pseudoscalar condensates and the two-flavor superconducting 
phase has also been examined j2Sl; although this remains to be applied to the charge 
neutral equation of state. 

The NJL model was originally developed to describe interactions between 
structureless nucleons j2^. It has been shown that this approach does not lead to nuclear 
stability (23 12S1 12Z| ■ The idea of the NJL model in its present form is to describe QCD 
at low energies by assuming point-like interactions between quarks. This approximation 
is justified by the fact that gluon degrees of freedom should be of relatively minor 
importance at low energies. Therefore it is sufficient to construct a model where the 
gluon interactions are "frozen in," meaning they are only present implicitly through the 
couplings of the model. The general form of the NJL Lagrangian density is given by, 

Cnjl = ^{i^ -m)^ + Y, G^{i,T^i,)\ (LI) 

a 

where "0 and if) are the quark fields, m is the current quark mass and Ga are the coupling 
constants associated with the various interaction channels. The interaction Lagrangian 
density can be expressed in a variety of equivalent forms using Fierz transformations 
PH] . Here we decompose it into qq and qq interaction terms as follows: 



^^11 = - G.(^75r^)' (1.2) 

-G^{,l^ri^f-G,{i,Yr^f (L3) 

+ a(^75C^r2/5^V^^)(V^^C-S5r2/3^V^) (L4) 

+ Ga(V^7MC'rr2/?V)(V^^C-S^r2r/?^V') (1-5) 



In order of appearance the above terms correspond to the scalar, pseudoscalar, vector- 
isoscalar and vector-isovector qq interactions, and the last two terms correspond to 
the scalar and the axial vector diquark channels. The color matrices are given by 
[3^ = a/3/2 A"^ {A = 2,5,7), indicating that these are attractive, color anti-symmetric 
diquark channels. In the NM phase the diquark interactions will lead to color singlet 
nucleons, and in the QM phase to color superconducting pairs. Indeed, QCD supports 
the idea that there is a strong attraction in the color anti-symmetric fiavor singlet scalar 
channel, leading to the formation of condensed scalar diquarks (the so called 2SC phase) 
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[U]. In our model these condensed pairs should arise from the same interaction (1.4) as 
the scalar diquarks in the nucleon constructed in the Faddeev approach. || 

Because the NJL model is non-renormalizable, the model is only fully specified 
after the choice of regularization scheme ^9]. In our work quark confinement is 
ensured mathematically through the introduction of a finite infra-red cut-off (Air) in 
the proper time regularization scheme pT| EH! • In this method, the unphysical quark 
decay thresholds are eliminated PP, as also found in the Dyson- Schwinger approach |3Uj . 
Of course, in the deconfined phase the infra-red cut-off will be set to zero |3T] . 

2. Nucleons and nuclear matter 

To describe NM we first consider the nucleon and its internal degrees of freedom. 
Incorporating the quark substructure of the nucleon through the Faddeev approach will 
allow us to examine how the nucleon properties change with density. The homogeneous 
Faddeev equation for the vertex function in the nucleon channel {J,T) = (1/2, 1/2) 
has the form Tn = ZUn Tn, where Z is the quark exchange kernel and Un the product 
of the quark and diquark propagators For our present, finite density calculations we 
restrict ourselves to the static approximation [32j, where the momentum dependence of 
the quark exchange kernel Z is neglected and a static parameter, c is instead introduced 
to reproduce the main features of the exact Faddeev calculation H^. Then Iliv becomes 
effectively the quark-diquark bubble graph given by 



where Tg and rj^'^ refer to the scalar and axial vector components of the diquark t- 

matrices respectively % , and S{k) is the constituent quark propagator. The nucleon 

mass follows from the requirement that the Faddeev kernel K = ZUn has eigenvalue 1. 

In a similar way the mass of the Delta resonance can be calculated from the pole position 

in the (J, T) = (3/2,3/2) channel, where only the axial vector diquark contributes. 

II Note that the SU{3)c color symmetry is broken in the 2SC phase - the pairing is between just two 
colors and the third color remains unpaired. We do not assume axial vector diquark condensation in 
QM, as this would break further symmetries. 

% It is possible to describe the nucleon with just the scalar diquark channel and this has been considered 
elsewhere [33] ■ As in Ref. we use the approximate 'constant + pole' forms of the diquark t-matrices. 




(2.1) 



(2.2) 
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Here we wish to use a more accurate description of the nucleon and delta masses, 
including also the pion cloud contributions This will enable us to consider a 

quantitative comparison of the scalar pairing found in the confined and deconfined 
phases. We take into consideration the self-energy contributions corresponding to the 
— i> A^vr — i> N and — Att — processes for the nucleon and the A Avr A 
and A Ntc A processes for the A [HS]- The corresponding self-energies are, 

= O-JlN + (^NA (2-3) 
Sa = O-AA + ^AAf' (2-4) 

where 

where ojbb' = {Mb — Mb') is the physical baryon mass splitting (e.g. una = 1232 — 939 
MeV), and iu(k) = \/k'^ ^ rn^ is the intermediate pion energy. For the ttBB' vertex we 
assume the phenomenological dipole form u{k) = A^/ (A^ -|- k"^)"^. The coefficients cbb' 
come from the standard SU(6) couplings (i.e. = 1, cna = 32/25 ISE])- The baryon 
masses including the pion loop contributions are then given by Mb = M^ + S^, where 
are the "bare" masses which follow from the quark-diquark equation. 

The size of the self-energy of the nucleon. Eat, is not precisely known. Calculations 
indicate it could be up to -400 MeV [13 EH EHl EHj • In the present work Sat is varied by 
changing the dipole cut-off. A, within physically acceptable limits, in order to investigate 
how the pion cloud of the nucleon influences the equation of state of the system. 

The form of the effective potential in the mean field approximation has been 
derived for symmetric NM in Ref.pT], starting from the quark Lagrangian, Eq. ()1.5|) . 
and using the hadronization method. For the present calculations the effective potential 
is extended to the isospin asymmetric case, because for neutron star matter charge 
neutrality and chemical equilibrium typically lead to an abundance of neutrons over 
protons. The effective potential can be written as follows: 

V _Kac + V^-^-^-Y^, (2.6) 

where 

f d'k e^M^ (M-ml (Mo - my 

= ''V w + ^2-^) 

+ The dipole cut-off A should not be confused with Ajr and Auv which are the cut-off parameters on 
the quark level in this model. 
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is the vacuum term. The Fermi motion of the nucleons moving in the scalar and vector 
mean fields gives rise to the term 

V)v = -2 5^ / 7^ - ^) (V^i + M%{M) - ^k^ + MUM)) (2.8) 

where Mi\f{M) is the nucleon mass in-medium, which is the sum of the bare mass and 
the pion cloud contribution. The relations between the chemical potentials, which are 
the variables of V^^, and the Fermi momenta appearing in Eq. ()2.8|) are given by 

/i„ = ^Mjj + kl^ + 3cuo - po (2.9) 

fip = ,jMj~hk\ + 3uo + po (2.10) 

The constituent quark mass, M, and the mean vector fields (c^o and po) in NM are 
determined by minimizing the effective potential for fixed chemical potentials *. The 
chemical potential of the (massless) electron in the last term of ()2.6p is fixed by the 
requirement of beta equilibrium as pe = Pn — Pp- 

For discussion of the phase structure of this model in Sect. 5, we introduce the 
chemical potentials associated with baryon number and isospin: 
1 1 

f^B = -{f^p + f^n); p/=-(pp-p„). (2.11) 

The parameters of the model are determined as follows. We choose Air = 285 MeV 

to be of the order of Aqcd- We calculate Auv, rn and so as to reproduce f-„ = 93 

MeV (through the matrix element for pion decay), = 140 MeV (through the Bethe 

Salpeter equation for the pion) and constituent quark mass at zero density, Mq = 400 

MeV (via the gap equation). The resulting equation of state is not very sensitive to the 

initial choices of Air and Mq. The parameter Gi^ is fixed to give the empirical binding 

energy per nucleon of symmetric NM {Eb = 17 MeV), and the parameter Gp is adjusted 

to the symmetry energy (04 = 32 MeV at po = 0.17fm~'^). The coupling in the axial 

vector diquark channel, Ga, is fixed by the mass of the Delta, since in this case the scalar 

diquark channel does not contribute. Then the coupling in the scalar diquark channel, 

Gs, is in turn determined by the nucleon mass. For the static parameter in the quark 

exchange kernel we use c = 450 MeV. This calculation of baryon masses is carried out 

for several initial choices of the dipole cut-off A, which controls the magnitude of the 

* In the actual calculation it is easier to minimize the energy density £ = V+Y^ ^J-aPa for fixed densities, 
after eliminating the mean vector fields as loq — 6Gui{pp + Pn), Po = 2Gpipp — Pn). 
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pion loop contributions to the nucleon mass, S^v. We will show the results for three 
different values of A, leading to S^v = -200 MeV, -300 MeV and -400 MeV. The 
resulting parameter sets are shown in Table 1. (The first set corresponds to the case 
where the nucleon and delta masses are reproduced without pion cloud contributions.) 



A 


Sat 


Sa 


















0.27 


0.30 


0.48 


0.70 


710 


-200 


-183 


0.21 


0.27 


0.30 


0.74 


803 


-300 


-266 


0.18 


0.24 


0.24 


0.76 


877 


-400 


-346 


0.15 


0.22 


0.19 


0.77 



Table 1. Parameters corresponding to different choices of Sat. We define — Gq,/G^ 
(a = a,s,u!,p). A, S^r and Sa are in MeV. 

The function Mn{M), which is needed in (12. 8j) to minimize the effective potential 
for finite density, is calculated by assuming that the ratio qa/ fn in ()2.5|) is independent 
of density. This is supported by a recent analysis of the pion-nucleus optical potential 
|4()j . which has shown that the pion decay constant in nuclear matter is reduced by 20%, 
which is the same as the quenching of qa derived from Gamow- Teller matrix elements 
|41j . The pion mass in the medium is constrained by chiral symmetry which leads to 
the relationjj m^(p)^ = ■ Mq/M. This small enhancement of the pion mass in the 
medium is taken into account in the calculation but its effect is not very important. 

The relationship between the nucleon mass and the constituent quark mass for the 
four cases of Table 1 is shown in Fig^ Note that in the region M > 200 MeV, which 
is most relevant for normal densities, both the slopes and the curvatures of the lines 
decrease in magnitude as the attraction from the pion loop increases. Thus the inclusion 
of pion loop effects on the nucleon mass leads to a reduced effective a NN coupling and 
to a reduced scalar polarizability in the medium |Tj. The self consistently calculated 
quark and nucleon masses are shown as functions of the density in FigO 

Concerning the pion exchange effects on the NM equation of state, we note that 
in addition to the term which can be incorporated into the nucleon mass there is 
also the familiar "Fock term" which originates from the Pauli principle. However this 
contribution is quite small when the short range correlations between nucleons in the 
spin-isospin channel are included ^21, and will be neglected here for simplicity. 
For the derivation see Ref.^j. The pion mass in this relation is defined at zero momentum. 
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Figure 1. The mass of the nucleon as a function of the scalar potential ($ — Mq~M). 
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Figure 2. The effective masses of quarks (lower lines) and nucleons (upper lines) in 
the hadronic phase for several choices of the nucleon self-energy. 

3. Quark matter 

In the most simple mean field description of quark matter, each energy level for the 
quarks is filled up to the Fermi energy. At the Fermi surface, only a small attraction 
between quarks leads to the formation of Cooper pairs, analagous to the phenomenon 
of BCS pairing of electrons in a superconductor. In QCD we anticipate that color 
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superconductivity will arise through the 3c channel, reducing the energy of the system 
through the condensation of color anti-symmetric pairs. The effective potential for QM 
thus has the form 

^^"^ = Kac + Vq + Va--^ (3.1) 

where the vacuum part, Kac, is the same as for NM (Eq. ()2.7|l ) except that the infra-red 
cut-off is zero in QM. Furthermore, 

— e(/i, - EQ{k)) {^^^ - EQ{k)) (3.2) 

a=u,a ^ 

describes the Fermi motion of quarks with chemical potentials /i„ and yU^, and Eqlk) = 



vMM-P^. The term ()3.2j) is analagous to Vn in NM, except that the quark mass, M, 
corresponds directly to the scalar field in the system, i.e., there is no scalar polarizability 
of the quarks. The term Va describes the effect of the pairing gap and is given by 

where e±{k) = a/ {Eqik) ± ^ib/^Y + A^, E± = \EQ{k) ± /iB/3|, and we introduced the 
chemical potentials for baryon number and isospin ff 
3 1 

fiB = ^ifJ-U + fJ-d); /i/ = -(/i„ - yUrf), (3.4) 

which corresponds to (j2.1H) in the NM phase. The electron chemical potential is 
determined from beta equilibrium as fie = f'd — fJ'u = —'^f^i- 

The gap A and the quark mass M are determined by minimizing the effective 
potential for fixed chemical potentials. Our results, discussed below, show that M is 
quite small in the QM phase, i.e., we have almost current quarks. 

Note that the vector-type interactions are set to zero in QM, even though in the 
description of NM the vector mean fields are clearly important. This assumption, which 
has been made implicitly in almost all investigations of high density quark matter, is 
supported by the discussions of vector meson poles in Ref. |31j. It is also supported by 
recent arguments related to the EMC effect 45 , which show that in the high energy 
region, where one has essentially current quarks (as in the present high density case), 
the mean vector field must indeed be set to zero. 

ff We mention that in principle one needs a further chemical potential for color neutrality (/is) in QM. 
However, for the 2-flavor case /is turns out to be very small |15l 1431 HI) . 
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Figure 3. Phase Diagrams with rs{NM) — rs{QM). The nucleon self-energies (Ejv) 
are 0, -200, -300 and -400 MeV respectively. Note that in each region the sign of the 
charge density is indicated. Globally charge neutral matter occurs at the boundary of 
the positively and negatively charged regions. 



In the following discussions, we distinguish the normal quark matter (NQM) phase, 
which is characterized by A = 0, from the color superconducting (SQM) phase (A > 0). 
Note that the value of Gg = rsG.,^ controls the outcome of this minimization through the 
last term in In these calculations we focus on the case where the pairing strength 

(r^) in QM takes the same value as in NM, namely the value required to obtain the 
correct nucleon mass after the pion cloud contribution is taken into account. 
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To construct a phase diagram in the plane of the chemical potentials, as in Fig. 3 which 
is discussed below, we compare the effective potentials for the NM, NQM and SQM 
phases at each point. The phase with the smallest effective potential (largest pressure) 
is the one which is physically realized. At each point we calculate the baryon and charge 
densities, 

1 dV dV 

where a = NM, NQM or SQM. The actual equation of state as a function of baryon 
density is then constructed according to global charge neutrality. If charge neutrality can 
be realized within one phase (for example the NM phase), one simply moves along the 
charge neutral line in the phase diagram (for example the line +NM/-NM in Fig 121) as 
the baryon density increases. When a phase transition occurs, it is necessary to construct 
a mixed phase, which is composed of positively and negatively charged components 
belonging to two different phases jlHI • A charge neutral mixture of NM and QM (where 
QM refers to either NQM or SQM), for example, is characterized by the volume fraction, 

y^A^ = f4 7) 

He Pc 

which ranges from to 1, as the density increases from the point of pure NM (where 
p^A^ = 0) to the one of pure QM (where p^*^ = 0). 

The baryon and energy densities for the mixed phase are then expressed by the 
volume fraction as follows, 

^ ^NMpNM ^ (1 _ ^NM^pQM 

^ ^NM^NM ^ _ ^NM^^QM (4 

Note that the components of the mixed phase have equal pressures {P^^ = = pQ^'^ 
at each point {fiB,Pi) on the phase boundary. In this way one moves along the phase 
boundaries between NM and QM while is increasing, until one comes to the point 
where charge neutral pure QM is realized {x^^^ = 0) 

In practice, our procedure is as follows. We first find the point where the effective 
potential for charge neutral NM becomes equal to the one for QM. At this point 
■y^NM _ -^^ since p^^^ = 0. From this point we incrementally increase either the 
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neutron or proton density (depending on whether the transition is in the direction 
of increasing neutron and/or proton density). For example, the transition from NM to 
SQM {ts = 0.27) is in the direction of increasing neutron density. For each neutron 
density, we determine the value for proton density required to ensure that we are on the 
phase boundary. Next we calculate the charge densities of each phase and the resulting 
volume fraction. The volume averaged properties of the mixed phase are then related 
to its component phases by Eqns. (j4.8p and (j4.9|) . This process continues until we 
encounter the point where the QM phase becomes charge neutral {pf^^ = 0), and thus 
the volume fraction x^^^ goes to 0. From here the remainder of the equation of state 
will be the pure charge neutral QM phase. 

5. Results 

The equations of state for this model exhibit phase transitions from the confined quark- 
diquark states employed in the description of NM to a phase with condensed quark pairs 
in the form of color superconducting QM. The key point in our present work is to equate 
the pairing strength for the color superconducting pairs in QM with the scalar diquark 
interactions inside the nucleon. Depending on the amount of attractive contributions of 
the pion cloud to the nucleon mass, this gives rise to a series of phase diagrams shown 
in Fig 121 corresponding to the four cases of Table 1. 

The four diagrams in FigOl are shown for increasing pion cloud contributions to 
the nucleon mass, i.e., decreasing scalar diquark interactions characterized by the ratio 
Tg. Starting from the first diagram in Fig|21 the regions occupied by SQM become 
smaller while those of NM and NQM phases become larger. In the first phase diagram, 
where the pion cloud is not included and the whole attraction within the nucleon is 
attributed to the diquark correlations, the NM phase is expelled almost completely and 
overwhelmed by the SQM phase even in the low density region. Since this is clearly 
unphysical, we can conclude that the naive quark-diquark picture of the nucleon leads 
to a phase structure which is in conflict with empirical evidence. Hence some attraction 
within the nucleon must be attributed to the pion cloud. 

The second diagram in Fig. 3, which corresponds to a mass shift of S^v = —200 MeV 
from the pion cloud, involves a mixed phase -NM/+SQM before the system undergoes 
a transition to the pure SQM phase. The phase boundary in this case is rather short. 
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Figure 4. Equations of state and neutron star masses for the choices of scalar 
attraction described in the text. 

indicating that the pressure is almost constant during the phase transition. In the 
third diagram (E^v = —300 MeV), the NM region extends to larger values of /i^, and 
the charge neutral phase boundary shrinks almost to a point. The fourth diagram 
(S^f = —400 MeV) involves a mixed phase +NM/-SQM before the pure SQM phase is 
reached. 

The charge neutral equations of state for these four cases are shown on the left hand 
side of FiglH In the low density region (0 > > 0.001 fm"'^) we use the equation of 
state of Negele and Vautherin jj^^ to describe the neutron star crust. On the NM side of 
FigEl the equation of state becomes softer with increasing pion cloud contributions, Sat. 
This can be understood from our discussions in Sect. 2. In particular, the inclusion of 
the pion cloud leads to a reduction of the effective aNN coupling in the medium, which 
must be balanced by a smaller vector coupling (the parameter in Table 1) in order to 
maintain the correct binding energy at the nuclear matter density 0.17 fm~^. Note that 
the actual saturation point for the model moves to somewhat higher densities (0.14, 
0.19, 0.21 and 0.25 fm'^) when S^v increases (0, -200,-300 and -400 MeV respectively). 
As noted in earlier works the only parameter in the present model that may be adjusted 
to give the saturation point is With this in mind these results are reasonably 
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close to the empirical value of po = 0-17 fm~'^. 

On the QM side all cases considered here give almost identical equations of state, 
indicating that the relationship between density and pressure in charge neutral SQM is 
not sensitive to the value of (or indeed to the value of A). However, the transition 
densities are sensitive to r^. In all cases considered here, the pressure variations in the 
mixed phase are rather small, indicating that our construction, based on the Gibbs 
criteria of phase equilibrium, gives similar results to a naive Maxwell construction 
between charge neutral NM and SQM. 

Through the Tolman Oppenheimer Volkoff (TOV) equations |1H] any equation of 
state specifies a unique set of non-rotating relativistic stars. The right hand side of 
FiglU illustrates the solutions to these equations for each of the equations of state on 
the left hand side of Fig0] We see from these figures that the inclusion of pion cloud 
contributions to the nucleon mass, leading to a decreasing strength of scalar diquark 
pairing, has significant effects on the properties of neutron stars. Because of the softening 
of the equation of state, the maximum masses for pure hadronic stars are reduced from 
2.4 solar masses to 2.0, 1.8 and 1.7 solar masses as S^r increases (and decreases 
accordingly) . 

The equations of state with phase transitions to SQM produce plateaus in the 
central density vs mass curves in FigHJ In the case of = 0.22 only the mixed phase 
can be reached inside a star, since a region of negative slope in this plot corresponds to 
unstable solutions to the TOV equation. However, stable hybrid stars with quark cores 
are possible in the case of Tg = 0.27 and = 0.24. In these configurations there are 
additional sets of stable solutions at higher densities, which in the literature are refered 
to as twins jlHI (since for these configurations there can be stars that have the same mass 
but different radii, as illustrated in Figl^I). Qualitatively similar results are also found 
in Ref . |50j . It has recently be shown that the second set of solutions are indeed stable 
and may give rise to an observable signature for the occurance of phase transitions in 
compact stars 

However, in the case of = 0.27 the masses are too small to allow for observed 
pulsar masses, which are typically about 1.4 solar masses. For Vg = 0.24 and = 0.22 
the maximum masses are approximately 1.4 and 1.5 solar masses, respectively. One 
possibility that may give rise to more massive hybrid stars in this model is to delay 
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Figure 5. Mass radius relationship of neutron stars and hybrid stars for the choices 
of scalar attraction described in the text. 



the phase transition to QM. If hyperons were included, for example, then the transition 
may be shifted to higher densities, since the equation of state in the hadronic phase 
would be softened independent of the value of the pairing strength r^. We note that 
the correlation between the transition densities and the maximum neutron star masses 
shown in FigE] follows the phenomenological discussions on phase transitions given in 
Ref.[SSl- 

It is interesting to note that these phase transitions to QM give rise to plateaus in 
the neutron star masses. This phenomenon may be the reason that so many observed 
neutron stars lie within such a narrow mass range |S2|. In addition. Fig. El shows that 
the radii of the stars may be reduced by 2 - 4 km if the pion cloud contribution to the 
nucleon mass increases (and the scalar pairing interaction decreases accordingly). 

6. Discussion 

We have used the NJL model, supplemented with a method of regularization which 
simulates confinement, to calculate consistently the properties of nuclear matter and 
quark matter. In doing so we have related the quark interactions in the confined phase, 
to the quark interactions in the deconfined phase, where color superconductivity is 
assumed to arise. Our principal finding is that the scalar pairing between quarks within 
the nucleon and in QM may be equated, if the attraction within the nucleon is attributed 
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not only to the diquark interactions but also to the pion cloud. Since the attraction in 
the scalar 3c channel between (almost) current quarks at high densities or energies can 
be derived directly from QCD, this result lends some support to the Faddeev approach 
to the nucleon, since the quark-diquark picture may be characterized by not only the 
same type of pairing interaction but also the same strength as we expect to find in 
the high density QM phase. This is an important feature of nucleon dynamics and is 
relevant to any finite density studies that incorporate nucleon structure. 

By including the pion cloud contributions to the nucleon mass, we found that the 
equation of state of nuclear matter becomes softer, reducing the neutron star and hybrid 
star masses significantly. 
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